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Semiconductor microcavities with artificial single-photon emitters have become one of the backbones of 
semiconductor quantum optics. In many cases, however, technical and physical issues limit the study of optical 
fields to incoherently excited systems. We analyze the model of an incoherently driven two-level system in a 
single-mode cavity. The specific structure of the applied master equation yields a recurrence relation for the 
steady-state values of correlations of the intracavity field and the emitter. We provide boundary conditions that 
permit a systematic solution which is numerically less demanding than standard methods. The method allows 
us to directly infer reasonable cutoff conditions from the system parameters. Different cavity systems from 
previous experiments are analyzed in terms of field correlation functions which can be measured via homodyne 
correlation measurements. We find that nonclassical correlations occur in systems of moderate quantum-dot- 
cavity coupling rather than strong coupling. Our boundary conditions also allow us to derive analytical results 
for the overall quantum state and its higher-order moments. We obtain very good approximations for the full 
quantum state of the field in terms of the characteristic functions.lt turns out that for every physically reasonable 
set of system parameters the state of the intracavity field is nonclassical. 

PACS numbers: 42.50.Pq, 42.50.Ct, 37.30.-i-i, 03.65.Fd 


I. INTRODUCTION 

The basic structure of a two-level system (TLS), located in 
a quasi-resonant single-mode cavity, shows remarkable quan¬ 
tum properties. Already in the regime of weak coupling be¬ 
tween the intracavity field and the TLS, the emission rate of 
the latter is increased by the so-called Purcell factor ITJ. For 
stronger coupling, the energy eigenvalues of the system are 
drastically changed, resulting in dressed states (polaritons) 
and Rabi splittings El, which are also known from strong 
atom-laser interaction Em. The mentioned phenomena are 
pure quantum effects indicating the quantum nature of the 
systems under study. Thus, such a setup also constitutes a 
versatile source for nonclassical light. Antibunching and sub- 
Poisson photon statistics were shown in experiments with ions 
in optical cavities SID. For details, see also El El. 

More recently, semiconductor microcavities became a fo¬ 
cus of research. In these systems, excitons in quantum dots 
adopt the role of the TLS @. Semiconductor microcavities 
are much more complex than their atomic counterparts, as 
both the cavity and the quantum dot are embedded in an inter¬ 
acting medium. In particular, the dissipation rates of dot and 
cavity are usually dominant, limiting the possibility of strong 
coupling. Nevertheless, antibunching and sub-Poisson pho¬ 
ton statistics could be demonstrated for quantum dots them¬ 
selves Eol HH as well as inside microcavities lfT2l . Fur¬ 
thermore, Rabi splitting was achieved in some realizations of 
semiconductor microcavities CMB]. 

Another peculiar aspect in semiconductor microcavities is 
the difficulty to drive them with coherent light. This is due 
to multiple factors such as intensive light scattering off the 
sample geometry, resonance frequencies unfavorable for cur¬ 
rent laser technologies, and others im. As a consequence. 
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incoherent excitation is a broadly applied method in semicon¬ 
ductor optics, e.g., via far-detuned photoluminescence ifT^ 
or electroluminescence M- Theoretical descriptions of the 
quantum fields emitted from such systems are of great rele¬ 
vance to semiconductor quantum optics. For excitons behav¬ 
ing like bosonic particles, it has been shown that the regimes 
of weak and strong coupling depend sensitively on the system 
parameters ini. 

An early treatment of the steady-state properties of inco¬ 
herently driven TLSs in cavities was given by Agarwal and 
Dutta Guppta im. They applied continued fraction methods 
to obtain solutions and compared them to full numerical cal¬ 
culations. Nonclassical phenomena of the radiation field had 
been studied on the basis of first- and second-order moments 
of the photon number statistic by inspection of Mandel’s Q 
parameter m. This is a first important step, which gives in¬ 
sight in the sub-Poisson photon statistics. Later on, more gen¬ 
eral criteria for the nonclassicality of light, such as matrices of 
moments of the photon number operator ESI, had been intro¬ 
duced. This method was further generalized to a full charac¬ 
terization of the quantum properties of light through higher- 
order moments of two noncommuting observables ll^ l22l . 
Alternatively, general nonclassicality tests can be based on 
the characteristic function of the radiation field ll^l24ll . Such 
methods have not been applied yet to the systems under study. 

In the present paper we study the quantum properties of an 
incoherently pumped TLS in a microcavity. Due to the inco¬ 
herent dynamics, only specific correlations couple with each 
other. This yields closed infinite sets of coupled equations. 
In the steady-state case, correlations between the quantum 
dot and the intracavity field follow from recurrence relations, 
which are numerically less demanding than previous meth¬ 
ods. We will solve these equations by applying boundary con¬ 
ditions following directly from the recurrence relations. The 
boundary conditions not only yield the necessary criteria for 
the solution but also an analytical proof of convergence. The 
structure of the solution allows us to determine an appropriate 
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cutoff for numerical calculations, which solely depends on the 
system parameters. The boundary conditions directly imply 
properties of the quantum state of light. In particular, the state 
of the intracavity field can never resemble a thermal one under 
realistic physical conditions. Based on general moments cri¬ 
teria for nonclassicality, we analyze various quantum effects 
of the system for experimental parameters of different cavities 
studied within the last decade. It turns out that the system with 
moderate quantum-do-cavity coupling shows stronger signa¬ 
tures of nonclassicality than for stronger coupling. Finally, we 
apply the boundary conditions to approximate the character¬ 
istic function of the quantum state of the intracavity field with 
controlled errors. In this way we prove the nonclassicality 
of the latter for any incoherently pumped quantum-dot-cavity 
system. 

The article is structured as follows. In Sec. we spec¬ 
ify the system under study and give the corresponding equa¬ 
tions of motion. Section|In]is used to introduce our ansatz for 
solving the steady-state case using recurrence techniques. We 
also derive the above mentioned boundary and initial condi¬ 
tions there. In Sec.|IV[ we test various nonclassicality criteria 
with parameters describing realistic systems applied in recent 
years. Section [V] deals with the characteristic function of the 
intracavity field, including a general treatment of the asymp¬ 
totic behavior of the recurrence relation. In Sec.|VI[ we give a 
summary and some conclusions. 


II. SYSTEM 


A scheme of the model under study with the considered 
processes is depicted in Fig. [T] The quantum dot in a micro¬ 
cavity is described by a TLS with transition frequency a;2i 
coupled to a mode of the electromagnetic held inside the cav¬ 
ity with frequency ojc = u!2i+6. Applying the dipole approxi¬ 
mation and the rotating wave approximation, and shifting into 
the frame rotating with a;2i, we obtain the Jaynes-Cummings 
Flamiltonian ll 25 l 

H = hS a^a + hg ^a^Ai2 + ^210^ . ( 1 ) 


Here, a and are the annihilation and creation operator of the 
held mode, respectively, Anm = l^)(w| (m^n = 1,2) are 
the atomic operators of the TLS, and g denotes the coupling 
strength. 

The dynamics of the system under incoherent pumping is 
governed by a master equation for the system density operator 
p El, which reads as 


df 


P = 


i 

n L 


HJ 




( 2 ) 


with pumping strength p, quantum-dot spontaneous-emission 
rate F, and cavity decay rate k. The relaxation 
processes are incorporated via the Lindblad operators 
Cj,{p) = [Xp,X^] + [X,pX^]. 

In the following, we write down all correlations in normal 
ordered form for the sake of clarity. The structure of Eq. 0 



FIG. 1. (Color online) Scheme of a quantum dot (two-level system) 
in a (semiconductor) microcavity. The left mirror is supposed to re¬ 
flect perfectly, whereas the right one is partially transparent giving 
rise to a loss of the cavity mode. 


has the effect that only those moments (^A2iA\2(A'^a'^i 
{k,l,m,n G N) that create or annihilate the same num¬ 
ber of excitations overall (that is, both intracavity field and 
TLS) couple to each other. Furthermore, due to the inco¬ 
herent nature of the pumping process, those moments with 
excess creation or annihilation vanish in the long-time limit 
as they do not couple to the only inhomogeneous correlation 
1 = (A21 ^ 520 ^ °a°). There remain three types of indepen¬ 
dent, nonvanishing moments, which we abbreviate as 

/„ = (at"a"), ( 3 a) 

= (i22at”a"), ( 3 b) 

= (i2iat"a"+i). ( 3 c) 


The corresponding set of coupled equations of motion can eas¬ 
ily be derived from (|^ and reads 


-/ = 
df " 

dt 

-R - 


—nuln — 2nglm[Rn-i], 

-a„Bn + pin + 25 Im[i?„], 

(iS 'Yn) Bn 

— ig[{n + l)Bn + 2 Bn+i — In+i]-, 


( 4 a) 

( 4 b) 


( 4 c) 


where we dehned 7„ = (F -f p -F K{2n + 1)) /2 and (T„ = 
F -F p -F UK. On a side note, we state that for certain limit¬ 
ing conditions, analytical solutions can be found; e.g., in the 
steady state for k — 0, we find a thermal state for the intra¬ 
cavity moments In the following, we consider all dissipa¬ 
tion rates and g to be finite but nonzero, and define this as a 
physically reasonable set of system parameters. 

It should be noted that, in general, semiconductor quantum 
dots also experience strong nonradiative dephasing. This can 
be easily included by adding a Lindblad term ^C^^^{p). It 
couples only to the moments yielding a different value for 
7„, which reads 


7n = 7n + 


( 5 ) 
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For our purpose, it does not change any of our qualitative 
statements and the amplitude of the considered correlations 
varies only slightly. Furthermore, for the realistic cavities con¬ 
sidered in Secs.|IV|and[V]no values for Fd were given. Hence, 
we discard these dephasing effects throughout the paper. 


For a„+i > 0, this inequality is trivial to fulfill. However, for 
a„+i < 0 , it yields a bound for the growth of the moments /„ 
as 


^n+l ^ l^n 
In Oiji+l 


(14) 


III. ANSATZ FOR STEADY STATE SOLUTIONS 

We are interested in the long-time solution (steady state) 
of the set of equations Q. In this case, Eq. Q reduces to 
a homogeneous and coupled set of algebraic equations. The 
specific structure of these equations makes it feasible to for¬ 
mulate a recurrence relation for the moments /„ of the form 


Now, for sufficiently large n we may only take into account 
the leading terms of the coefficients, which read (see Ap¬ 
pendix]^ for details) 


Q^n +1 

/3„ 




2k 


In-\-2 — ^n+lIn+1 “f jlnlnt 

with the recurrence coefficients a„, /3„ given by 


( 6 ) 


After inserting these into Eq. ([T4J, we obtain 


In+l ^ Ig P 1 _ ^ 1 

7- — Si 

In n n 


(15) 

(16) 

(17) 


— 


/3n = 


Un ( 2p riK 

1 + A„_i \ 

(7) 

with the positive constant ^ = 2g^p/ kP . As this ratio goes to 

2,K (7n—l 

A„_1 )' 

zero, we can set as a second boundary condition 


{n + l)pan+i 


( 8 ) 

lim In = 0 . 

(18) 

2 k On ’ 


n—>oo 


where we set A„ = ( 2 p^ 7 „) /k (S'^ + 7 ^). 

In order to solve this recurrence relation, we require two 
independent boundary conditions. Once found, relation 
yields the full steady-state solutions for all moments Bn, 
and Rn- The first condition, related to the normalization or 
completeness relation, reads 


lo = (i) = 1 . 


(9) 


The other one, usually applied by truncating the endless hi¬ 
erarchy of coupled equations, would follow from the bound¬ 
edness of the density operator g. As lim {n\g\n) = 0, we 

n—FOO 

may set (nolpl’^o) = 0 , thus cutting off the coupling between 
the different moments for n < uq and n > uq, yielding a 
closed system of equations. The error from this approach can 
be adjusted by increasing uq- The resulting boundary condi¬ 
tion reads as 


(no|p|no) =— (: exp(- 

no! 

- (- 1 )/ - 
tt-q! k\ 


, , , ,1 Ino+k — 0 , 

tt-q! k\ 


■a^a) :) (10) 

( 11 ) 


where (: • :) denotes normal ordering. In general, this crite¬ 
rion is difficult to apply. Another possible formulation of the 
second condition is found using the fact that for any order 
n, the moment In is positive semidefinite. Therefore, using 
Eq. (j^, we find 


0 ^ Oin+lln+l “f Pnin- 

Since Pn > 0 holds for all n, this is equivalent to 

In. . Q^n-t-1 


> 


^n +1 


Pn 


( 12 ) 


(13) 


A few things should be noted. Eirst, Eq. ( [TSl l is by no 
means a general result for quantum states of bosonic systems. 
Quite the opposite, the moments for both a coherent state 
(/“h = If) with Ji > 1 and a thermal state (I*"™ = n!/f) 
with /i > 0 diverge. In particular, the latter is relevant as 
it would be the state in either of the (not physically reason¬ 
able) cases of At —0 or p —00 (see Appendix |B|). The other 
way around, we directly conclude that neither of these states 
is reached in this setup. Second, in the Appendix |A| we also 
estimated a lower bound of the ratio In+i/In, approaching 
again ^/n for large n. Hence, the right-hand side of Eq. (17 1 


is not just an upper bound, but a very good approximation for 
the large-n behavior. We will use this fact when we consider 
the characteristic function of the intracavity field. Third, as a 
consequence of our second note, we observe that neither 5 nor 
F play a role in the asymptotic behavior. 

The boundary condition ( fTS] ) can be utilized to calculate 
a more practicable condition for Ii. Iteratively applying the 
identity (j^, we obtain 

In+2 = Cn+2ll + Dn+2lo- (19) 

The coefficients Cn +2 Dn +2 obey recurrence relations 
like (j^, namely 

Cn+2 — ^n+lCn+1 “t“ PnCn, 

Rn+2 — ^n+lRn+1 “f PnRn, 

with initial conditions 


Co = 0, 

Do = l, 


Ci = l, 

£>i = 0 . 


( 20 ) 


Applying the conditions Iq = 1 and lim In = 0, we find 

n—>-00 


/i = lim — 


Dn 

Cn 


( 21 ) 
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The calculation of the moments Bn and Rn is straightforward 
using the steady-state form of Eq. It should be kept in 
mind that for numerical calculations, a truncation has to be 
done in Eq. ( |2T] ), where we have to set /jv = 0 for an ap¬ 
propriate order N and only moments up to this order can be 
calculated from Eq. (|^. 

Let us consider the numerical complexity of our method. 
Since the initial conditions ( |20| are known and the and 
/3„ are analytical functions in n, the recurrence relations for 
the Cn and are of the order of 0{N). Using the recur¬ 
rence relation Q for the /jv together with the calculated Ii, 
our approach yields the same complexity of 0{N). Let us 
compare this to the standard approach for solving the steady 
state problem; cf. Il26l l27ll . The main step involves the in¬ 
version of the matrix of coefficients following directly from 
the master equation (|^ governing the dynamic of the density 
matrix. Truncation at N means a maximum of N photons in 
the cavity, and thus 2{N -f 1) possible states. Eor the den¬ 
sity matrix, this yields [2{N + 1)]^ c>c iV^ rows. A matrix 
inversion based on the standard Gauss-Jordan elimination al¬ 
gorithm scales in computational effort with 0{M^), M being 
the number of rows. Thus we obtain for a photon number cut¬ 
off of N an 0{N^) complexity scaling. In contrast to this, our 
method reduces the complexity to 0{N). 

We can also use the explicit formulas ( [T4l i and ( [T7] i as well 
as the lower bound from the appendix, Eq. to infer an 
appropriate cutoff for the numerical calculations together with 
the corresponding error bars. Eirst, one derives the bounds for 
n, for which the dependences in Eqs. ( [T5| ), ( [T6| ) become dom¬ 
inant. Then we use the upper and lower bounds to determine 
the value for n, for which In+i/In is within the desired e 
neighborhood. Applying Eq. to Eq. ( |A5[ ), with the as¬ 
sumption of Cn+i > 0, we find 


Rji+i 
Cn+l — 


Al 


£-a„+i 



e-a„+i 


D, 

C 


<Ii< 


Dn+l — 

Cn+l — 


4r. 

-«n + l 
— dn+l 


D. 

C 


. ( 22 ) 


Comparing with Eq. ( |2T] l and the approximation in Eqs. ( [T5| l 
and ([T6| l, we see that the correction to the applied formula 
Eq. pi[ ) is of the order 0{n~^). 

As an example, we plotted in Eig. [^the first moment Ii, 
which is the intracavity field intensity, against the normalized 
coupling strength g/n. The parameters are the same as in the 
previous work of Agarwal and Dutta Gupta (cf. ifTSl ) and our 
method yields the same result. As the coupling increases, so 
does the intensity, and Ii ultimately reaches a saturated state 
for large coupling. The saturation values appear to be linearly 
dependent on the pumping strength. 

Another example for the use of the moments is the calcu¬ 
lation of the nth-order normalized intensity correlation of the 
intracavity field for zero time delav ll28l 


5(")(0) = 


In 

ir 


(23) 


Those normalized correlations are, in principle, measurable 
by higher-order generalization of the antibunching experi¬ 
ment by Kimble et al. Il29l , e.g., by time-multiplexing tech¬ 
niques IMI . and allow a comparison to well known states. Eor 
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FIG. 2. 7i = (cRa) for resonance 5 = 0 as a function of gj k. 
The normalized quantum-dot decay is T/k = 1 and the normalized 
pumping strength is varied from p/k = 0.5 to 2 in steps of 0.5, 
where a longer dashing corresponds to a stronger pumping. 



example a coherent state and a thermal state, as mentioned 
above, have moments = If and = n!/". Hence, 

their nth-order normalized correlations are equal to unity and 
n!, respectively. The intracavity field transits to a thermal state 
for p —> 00 , as can be seen from Eig. where we plotted 
/„/(/" n!) for n = 2,3,4 against the normalized pumping 
strength. We chose two parameter sets with normalized cou¬ 
pling gjn — 1 (blue/darker) and p/k = 5 (orange/lighter), 
each with T/k — 1 and i5 = 0. The curves tend for both 
sets and all orders to unity for increasing p/n. It should be 
kept in mind that albeit the system converges formally into a 
thermal state for p —> oo, it never reaches it. The transition 
is only possible since in Eq. ([T^ the upper bound ^ diverges, 
allowing the moments /„ to diverge. A rigorous proof of the 
thermal state limit is given in the Appendix [B] 



FIG. 3. (Color online) Renormalized correlation functions 
In/{Iin\) for parameters F/ft = 1, 5 = 0, and gj^, = 1 
(blue/darker) or p/k = 5 (orange/lighter). Depicted are the orders 
n = 2, 3, 4 with increasing dashing size, respectively. 
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IV. MOMENT BASED NONCLASSICALITY CRITERIA 

We proceed with an investigation of nonclassicality based 
on the moments extracted from the steady-state. Nonclassical 
light fields are quantum fields with properties which cannot 
be described by the classic Maxwell equations and are thus 
interesting for many applications. Agarwal and Dutta Gupta 
applied their method llTSl to analyze the Mandel Q parame¬ 
ter IIT 9 I of the intracavity field, 

( 24 ) 

(n) Ji 

Determining this parameter in experiments requires a mea¬ 
surement of the photon statistics as the amplitude of mean 
value and variance of the photon number need to be com¬ 
pared. In the following, we will focus on the more general 
correlation conditions according to ll20H22ll . There, we only 
compare moments with the same powers of field operators. 
Hence, their verification in experiments can be directly con¬ 
cluded from field correlations measured in homodyne correla¬ 
tion EH setups. 

The following moment based criteria are calculated for two 
exemplary microcavity systems, namely, micropillars and mi¬ 
crodisks. Those were realized in the last decade and give a 
good account of the realizable parameter ranges. The used 
parameters are taken from Khitrova et al. m and the two 
assigned parameter sets include the decay rates T and k of 
the quantum dot (TLS) and cavity, respectively, as well as the 
coupling strength g. Note that according to the used Hamilto¬ 
nian Q, all parameters are given as angular frequencies. The 
first set for the micropillars consists of T = 113 x lO^s-i, 
K = 276 X 10® s“^, and g = 122 x 10® s“^ and is denoted 
as set A. The second set, set B, for the microdisks is T = 
427 X 10® s-i, AC = 213 X 10® s-\ and p = 616 x 10® s"!. 
The remaining parameters, the pumping strength p and de¬ 
tuning 6 , can then be subject to the experiment performed 
on those microcavities, thus being the independent variables. 
However, we will set (5 = 0 since for our analysis the detuning 
is of less interest and we only vary p. 

Experimentally accessible criteria for nonclassicality can 
be formulated by considering principle minors of the matrix 
of moments for suitable basic operators; cf. Il20lj^ . In our 
case, we define 

/= (25) 

kGN 

fe,AG{0.1} 

Hence, the resulting matrix of moments is solely based on the 
moments defined in Eq. Q. Then the necessary and sufficient 
condition for a nonclassical state is that at least one principle 
minor of (: f^f :) becomes negative. However, we will only 
analyze certain minors related to particular nonclassical ef¬ 
fects, and thus the resulting conditions are just sufficient. The 
first condition deals with the intracavity field only and reads 

^ < 1. (26) 


The special case of n = 1 gives the well-known second-order 
correlation function for zero time delay, (0), and indicates 
sub-Poissonian photon statistics if = I 2 /I 1 < 1- 

InEig.g we plotted the condition for orders n = 1, 3, 5 
against the pumping strength, in each case for both parameter 
sets A and B. Eor all shown orders n, set A fulfills the condi- 



FIG. 4. (Color online) Condition I 2 n/In plotted for orders n = 
1, 3, 5 and parameter sets A (blue/darker) and B (orange/lighter) for 
varied pumping strength p. A longer dashing corresponds to higher 
orders n. Note that the order of magnitude of p is 10^ times higher 
compared to the other parameters; see beginning of Sec. 0 


tion ( |26l l below certain pumping strengths. Hence, the system 
undergoes a transition from the nonclassical to the classical 
regime for increasing pumping strength. Keeping the special 
case n = 1 in mind, this would be the transition from a sub- to 
super-Poissonian photon statistics. The set B does not fulfill 
the condition ( |2^ for any shown order n. Note that this does 
not rule out any nonclassicalities of higher orders. In fact, 
set B does fulfill the condition|26]for n = 10. 

The second condition reads 


B2n 

1 ^ 


1, 


(27) 


which is formally akin to condition ( |2^ but takes into account 
the intensity of both the TLS and the intracavity field. Con¬ 
dition <1^ is plotted against the pumping strength p in Eig. 
again for the sets A (blue/dark) and B (orange/light) and or¬ 
ders n = 1, 3, 5. Eor n = 1, the condition is not fulfilled for 
both sets A and B, but higher-order criteria are fulfilled for 
set A, again with consecutively lower values and for greater 
ranges of the pumping strength. This case again shows that 
the derived nonclassicality criteria are only sufficient and an 
unfulfilled criteria for some order n does not rule out nonclas¬ 
sical effects of higher orders. 

The third and last criterion examined here is 


Si/|i?o|"<l. (28) 


While the more general relation 

7?2n+l/ \Rn\^ < 1 


(29) 
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FIG. 5. (Color online) B 2 n /plotted for n = 1, 3, 5 and param¬ 
eter sets A (blue/darker) and B (orange/lighter) against the pumping 
strength p. A longer dashing corresponds to higher orders n. Note 
that the order of magnitude of p is 10^ times higher compared to the 
other parameters. 



FIG. 6. (Color online) Plot of Bi — |i?oP for set A (blue/darker) 
and set B (orangeAighter). Negative values occur for entanglement 
between the intracavity field and the TLS. 

indicates nonclassicality as the other criteria, Eq. ( [28] l is also 
a sufficient criterion for entanglement between the intracavity 
field and the TLS |[33l. In Fig. we plotted condition ( |28l l 
for our two example systems. Similar to the case of condi¬ 
tion ( |26l l, set B does not show entanglement, whereas set A 
fulfills the condition for certain pumping strengths. However, 
note the scaling of both the abscissa and ordinate. Negative 
values in Eq. ( [28] l only occur for very low pumping strengths 
and also attain only small values. 

The results so far indicate that set A shows nonclassical be¬ 
havior of different kinds, which reveal themselves via lower- 
order moment-based criteria. For set B, on the other hand, 
those same order criteria do not indicate a nonclassical charac¬ 
ter. In particular, it should be noted that an accurate measure¬ 
ment of higher-order field moments, such as Iiq as shown in 
Fig .[^ is currently out of technological reach (301 . This result 
is rather surprising as set B has a better quantum-dot-cavity 


coupling relative to the dissipation rates. We conclude that the 
creation of nonclassical light in incoherently driven quantum- 
dot-cavity systems may be favored by not too strong coupling. 
In semiconductor quantum optics, strong coupling, and espe¬ 
cially Rabi splitting, is considered a clear indicator of quan¬ 
tum light; see, in particular the argumentations in iflTHTSl . 
Our results, on the other hand, clearly show the occurrence 
of strong signatures of nonclassicality for moderate coupling, 
making this scenario favorable for applications. 

Detecting nonclassicality can be cumbersome, as there is an 
infinite hierarchy of sufficient but not necessary nonclassical¬ 
ity conditions (221. For some system parameters, lower-order 
moments do not reveal any nonclassical effect. Hence, in the 
following, we will investigate the characteristic function of the 
intracavity field, which includes the full information also on 
higher-order moments. It is important that the characteristic 
function can be directly sampled from experimental data (341 . 
Moreover, the characteristic function may uncover quantum 
effects more directly then moments criteria (351 . 

V. CHARACTERISTIC FUNCTION 

The characteristic function is one possible representation of 
the quantum state of light, thus containing all information on 
the state. The same information is contained in all the mo¬ 
ments In and, in this section, we derive their relation to the 
characteristic function of the intracavity field. Using this re¬ 
sult and the large-n behavior of the the asymptotic behav¬ 
ior of the characteristic function is found. With this we prove 
that the intracavity field is nonclassical for any physically rea¬ 
sonable set of system parameters. 

The characteristic function is defined as (^[^iJTll 

$(a) = D{a) , a e C (30) 

with the displacement operator Z?(q;) = exp(aa^ — a*a). The 
representation of $(a) in terms of the /„ is derived by eval¬ 
uating Eq. ( |30| , while bearing in mind that the off-diagonal 
elements vanish in the long-time limit due to the 

incoherent nature of the pumping. The result reads 

°° (— 1 )” 

= (31) 

n=0 ^ '' 

and is therefore phase independent, $(a) = <l>(|a|). For an 
actual calculation, we truncate the infinite sum at n = W in 
accordance to the truncation for the moments In = 0. 

The characteristic function is of special interest since it con¬ 
tains the whole information about the underlying state and ap¬ 
propriate nonclassicality criteria can be formulated (241 . The 
simplest criterion reads (231 

|^(a)| > 1- (32) 

In order to test Eq. ( |3^ , we plotted $(a) in Fig. |7]for the 
parameter sets A and B and a pumping strength of p = 
1 X 10^^ s“^. In the plotted range of |a| only set A fulfills the 









7 




FIG. 7. (Color online) Characteristic function for the parameter 
sets A (blue/darker) and B (orange/lighter) and pumping strength 
p = 1 X 10“ s-^ 


condition ( |3^ and therefore reveals nonclassicality. However, 
this analysis has its limits since only a restricted range of |a| 
can be displayed due to the truncation of the infinite sum in 
Eq. ® at n = N. Beyond these limitations, we can, how¬ 
ever, analyze the general behavior of the characteristic func¬ 
tion when using the asymptotic behavior of the moments 
The detailed derivation can be found in Appendix which 
gives the asymptotic relation 

/„+! = n>N, (33) 

n 

with ^ = 2pg^/ k^. Here, the order N is dehned as N-. = 
max(ni,712,ns), where ni, 77 - 2 , and 773 are the orders for 
which the approximations — a„+i oc 77 ^, /?„ oc n, and 
/„+! = In Pn/i—C(n+i) hold, respectively; cf. Eqs. ( |A3| l- 

It should be kept in mind that N can become consider¬ 
ably large and for actual calculations like in Sec. the re¬ 
lation ([ 3 ^ is not reasonable. Eor a visualization, we plotted 
in Eig.^lthe upper and lower bounds for the moments accord¬ 
ing to Eq. ( |A5| ) in Appendix [A| against the order 77 . We also 
added the asymptotic relation /?„/(—ctn+i) = ^/n as used in 
Eq. ( |3^ and the case were the next-to-leading order is kept 
for the coefficients a„+i and /?„; cf. Eqs. ( |A3| l,( |A4| ). The pa¬ 
rameters were chosen from set B with a pumping strength of 
p = 1 X 10^^s“^. Apparently, the asymptotic relation holds 
true for orders n « 10^, whereas the lower and upper bounds 
are already equal for n ~ 10. Hence, the moments /„ can al¬ 
ready be approximated by the upper bound when the asymp¬ 
totic relation still differs by orders of magnitude. 

While the above considerations are important for actual cal¬ 
culations involving any of the mentioned approximations, we 
can neglect them for the following theoretical analysis and 
solely rely on the fact that there is some N which fuffills the 
requirements of the asymptotic relation ( |3^ . Under these as¬ 
sumptions, we are now able to avoid the truncation but instead 
split the infinite sum in Eq. ( [3T] ) into two parts. In the hrst 
sum, we use the exact moments up to the order N, and in the 
second, we insert the asymptotic relation (|3^. The latter sum 


FIG. 8. (Color online) Plot of the lower (blue/darker) and upper 
(orange/lighter) boundaries for the moments. Also seen are the ap¬ 
proximated ratios /?„/(—On+i) were the linear order is kept (green 
dashed line) and the asymptotic relations hold (red dot-dashed line). 
The parameters were taken from set B and the pumping strength is 
p = 1 X 10^^ s“^. 


can be simplihed by iterative application of Eq. (H), leading 
to the following expression for <I>(a) 


N 


$(a) = ^ 


(-ir 

( 77 !)^ 




n=0 

iNiN-iy. 


e 


N 


t yffiop'-r. 


n—N-\-l 


( 771)3 


(34) 


With this expansion at hand, it can be shown that the asymp¬ 
totic behavior of the characteristic function is sufficient to 
prove nonclassicality according to the criterion ( [3^ for all 
physically reasonable sets of system parameters. At hrst we 
note that in general we do not know the coefficients of the hrst 
part of the expansion for the characteristic function as they 
follow from the full calculation. Hence, it is necessary for 
our aim to show that the asymptotic term increases faster than 
or any polynomial. Denoting x = |ap^, the asymptotic 
behavior can be reduced to 

(—x)" 

^’asymp(x) = V]\[ ^ ^ Tl / |\3 ) (35) 

n=N+l 


with Vn the positive prefactor. Note that x > 0 allows to take 
positive roots of x. As we are only interested in the order of 
asymptotic behavior, the values of both the scaling Vn and N 
itself are irrelevant. The latter would only imply a subtraction 
of a polynomial of order x^. We therefore set Vn = 1 and 
W = 0 to simplify our series to 

OO / s ^ 

‘^asymp(x) = ^ (36) 


Using Stirling’s approximation, we obtain 

( 77 !)^ « ^ ft! 27r773“^^""''^/^^(3r7)! (37) 
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and thus, after defining x = 


E 


n 


{-xY 

{n\Y 


V3 {-3xY^ 

2tt (3n)! 

n—1 ^ ' 


(38) 


The sum on the right hand side of Eq. ( |38| ) can easily be ob¬ 
tained by a sum of three scaled exponential functions with 
complex exponents, yielding 


^ (-3f)3» 
1 (3n)! 

n—1 ^ ' 



The addition of — 1 stems from the limitation of the series to 
n > 1. For a; ^ 1, both that constant term as well as the term 


exp(—3i) can be neglected. Inserting this result into Eq. (36 1 , 
we finally obtain 






1 


^37 


exp 


ia;3 


cos ( ^^a;3 ) (40) 


This is an alternating, diverging function with an envelope 

function going with exp (| a; 3), which will overcome any 
polynomial increase or decrease from the first part of the char¬ 


acteristic function in Eq. (34 1 or the approximation when set¬ 
ting = 0 in Eq. ( [36l l. Hence, we can state as our final result 
that the steady-state intracavity field in our system is always 
nonclassical. 


VI. SUMMARY 


We introduced an alternative technique for solving the 
steady-state problem of an incoherently pumped two-level 
system coupled to a single-mode cavity by using recurrence 
relations with appropriate boundary conditions. Due to the 
incoherent nature of the pumping, only three types of nonvan¬ 
ishing moments are needed, whose coupled equations of mo¬ 
tion lead to the recurrence relation in the steady state. Bound¬ 
ary conditions have been derived and an analytical proof of 
convergence of the moments has been given. Furthermore, 
we provided a direct estimate of the errors in dependence on 
the cutoff procedure to be used in numerical calculations. To¬ 
gether with the method of solution via recurrence relations a 
significant reduction of numerical effort was achieved. 

Nonclassicality criteria based on moments have been ana¬ 
lyzed for system parameters describing realistic microcavity 
structures. Nonclassical effects, such as sub-Poisson photon 
statistics for the intracavity field or its entanglement with the 
two-level system, occurred for certain parameters. The favor¬ 
able scenario for demonstrating quantum phenomena is that 
of moderate quantum-dot-cavity coupling and weak pump¬ 
ing. This is contrary to the often discussed claim that strong 
coupling is essential for revealing nonclassical signatures. In 
the limit of extremely strong pumping, the quantum state of 
the cavity field approaches a thermal one. For more realistic 
conditions, however, the state is far from showing a thermal 
statistics. 


We also examined the characteristic function of the intra¬ 
cavity field, which contains the full information about the 
quantum state of light, based on the asymptotic relations de¬ 
scribing the behavior of moments up to higher orders. This ap¬ 
proach renders it possible to show that the intracavity field is 
nonclassical for any realistic set of parameters. In conclusion, 
the techniques developed in our paper provide helpful tools for 
describing a two-level system in a cavity in the steady state. 
They offer a variety of possibilities for characterizing nonclas¬ 
sicality of the intracavity field and its entanglement with the 
degrees of freedom of the two-level system. 
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Appendix A: Lower bound and asymptotic behavior 

In addition to the upper bound for the ratio in 

Eq. ^ we find a lower bound being valid from a certain 
order on. The recurrence relation (|^ together with Eq. ( [I4| ) 
give 


^n+2 — ^n+l^n+1 T Yn^n ^ In+1 


Yn+1 


— C^n+2 

Hence, if there is a constant e > /?„+i/(—a„+ 2 ) 
have 


(Al) 
> 0, we 


^n+2 — ^n+l^n+l T Yn^n ^ ^^n+1- (A2) 


From the asymptotic behavior of the coefficients ttn+i, Pn, 
which follow from their definitions (|^an d|8), and read 


K 2 

/Sn ~ 


1 I ^(5 I 3r±P^ 

4g2 '' 2 ^ 2 K ) 


(A3) 

(A4) 


it follows that such a constant e does not only exist but that 
for every £ > 0 there is an order N from which Eq. ( |A2| ) 
holds. Note that for the special case of e < 1, the mo¬ 
ments In converge monotonically. N may be extracted from 
£ > 0n+il{—o-n+ 2 ) o r m ore e asily by using only the lead¬ 


ing terms in Eqs. (A3 and A4i and ^ = 2pg ^/as in the 


result (17 1 , giving N >^/e. 

Now, we can also formulate a lower bound and together 


Pn 


< ^ < 


Pn 


^ ^n+1 


^n+1 


n>i 


(A5) 


For sufficiently large n, the relation between /„ and In+i is 
well approximated by ^/n. 
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Appendix B: Asymptotic behavior for large pumping 


In case of very large pumping strengths p ^ T,K,g, 6, we 
can simplify the coefficients an+i and /3„ to 


ttn+l 


Pn 


P 

2k 4^2 
1 + n p 

2 K 


n + 1 
2 


(Bl) 

(B2) 


Applying again our upper-bound scenario from Eq. ( [T3] l, we 
obtain 


In _ 1 I^P 

In+I ~P ^ + 1 2p2(n-|-l) 


up 

2p2(n-h 1)' 


(B3) 


The latter approximation follows again for sufficiently large 
pumping. As in the limit p oo, the right-hand side can only 


be fulfilled for In+i 0, we can approximate the behavior, 
similar to the case above, by 

/„+!« +!)/„. (B4) 

\kp J 

Repeatedly applying this formula and taking into account that 
forp —oo, the value of n from which point on the asymptotic 
behavior is valid decreases, we find the relation of moments 
for a thermal state: 


In 


nl 



(B5) 


Consequently, the intensity Ii is then given by In the 
strict limit of p —>■ oo the moments = 0, Vn S N, and one 
obtains the vacuum state. 
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